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Abstract 

Nonlinear Schrodinger equation, complemented by a confining potential, pos- 
sesses a discrete set of stationary solutions. These are called coherent modes, 
since the nonlinear Schrodinger equation describes coherent states. Such modes 
are also named topological because the solutions corresponding to different spec- 
tral levels have principally different spatial dependences. The theory of resonant 
excitation of these topological coherent modes is presented. The method of mul- 
tiscale averaging is employed in deriving the evolution equations for resonant 
guiding centers. A rigorous qualitative analysis for these nonlinear differential 
equations is given. Temporal behaviour of fractional populations is illustrated 
by numerical solutions. 
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1 Introduction 



Nonlinear Schrodinger equation has recently attracted a great deal of interest, since it 
provides an adequate description for collective quantum states of trapped Bose atoms 
(see reviews [1-3]). This equation, as applied to Bose systems, is often termed the 
Gross-Pitaevskii equation, who first considered that physical application [4,5]. 

For the purpose of accuracy in terminology, it is worth mentioning the following. 
As is possible to show [6], the nonlinear Schrodinger equation is an exact equation for 
coherent states. While the correct meaning of the Gross-Pitaevskii equation [2] is that 
of an approximate mean-field equation for the broken symmetry order parameter. 

Coherent states of trapped atoms are the solutions of the nonlinear Schrodinger 
equation with a confining potential. Stationary solutions of this equation form a dis- 
crete set. Wave functions, corresponding to different energy levels, are called [6] coher- 
ent modes. Different coherent modes possess qualitatively different spatial behaviour, 
because of which they may be named topological coherent modes. Being the solutions 
of the nonlinear Schrodinger equation, these modes have nothing to do with collective 
excitations described by the linear Bogolubov - De Gennes equations. 

Bose-Einstein condensation of trapped atoms can be understood [3] as the conden- 
sation to the ground coherent state, that is, to the state with the lowest single-particle 
energy. In an equilibrium system, atoms always condense to the ground state. But, 
if an additional alternating field is switched on, with a frequency being in resonance 
with the transition frequency between two coherent energy levels, then higher topolog- 
ical modes can be excited, thus creating nonground-state condensates. The feasibility 
of such modes was first advanced in Ref. [7]. The alternating resonant field can be 
produced by modulating the confining potential of the trap or by imposing additional 
external fields. For instance, a rotating field for exciting vortices, can be created by 
multiple laser beams [8-10]. The resonant excitation can lead to the formation of sev- 
eral vortices [11] and, possibly, skyrmions in a spinor condensate [12]. The topological 
coherent modes have also been studied in Refs. [13,14]. Excitation of a dipole mode in 
a two-component condensate was observed [15]. The possibility of creating dark soliton 
states by means of the resonant excitation was studied [11]. As is found, dark soliton 
states are unstable with respect to their decomposition into several vortices. The states 
of multiple basic vortices are also more stable than a state of a single vortex with a high 
winding number [3,16,17]. Investigating the temporal behaviour [18,19] and collective 
excitations [20] of the formed coherent modes, anomalous dynamical phenomena were 
found [18,19] reminding a kind of critical phenomena. The origin of the latter has not 
been completely understood. 

The aim of the present paper is to develop the theory of the resonant formation 
of topological coherent modes and to provide an explanation of the dynamic critical 
phenomena. For this purpose, we give a rigorous stability analysis for the nonlinear 
evolution equations of guiding centers. We demonstrate that dramatic changes in the 
dynamic properties occur when crossing a saddle separatrix. The condition for the 
separatrix touching the boundary defines a critical line on the parametric manifold. 
This analysis is illustrated by the numerical solutions explicitly displaying a dramatic 
qualitative change in the dynamics of fractional populations when crossing the critical 
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line. 



2 Resonant Excitation 

Consider a system of Bose atoms, which have experienced Bose-Einstein condensation 
under conditions typical of experiments with trapped atomic gases [1-3]. The coher- 
ent state of such atoms is described [6] by a coherent field f satisfying a nonlinear 
Schrodinger equation with the nonlinear Hamiltonian 

H{ V ) = -^-V 2 + U{v) + NA\ V \\ (1) 
ZrriQ 

where m is atomic mass, U(r) is a confining potential, N is the number of trapped 
atoms, and A = Anh 2 'a s / 'm is the interaction strength, with a s being an s-wave scat- 
tering length. 

Topological coherent modes are defined [7] as stationary solutions ip n = y? n (r) of 
the nonlinear Schrodinger equation 

H(<Pn)<Pn = KVn • (2) 

Here n is a multi-index labelling the quantum states of coherent modes, and E n is 
a single-particle energy of the given coherent mode, which is normalized to unity as 
(ip n , (p n ) = 1. The excistence of the confining potential U (r) assumes that the spectrum 
{E n } is discrete, hence the set {y?„} of eigenfunctions is countable. 

It is worth mentioning that the solutions to a nonlinear Schrodinger equation do not 
necessarily form a basis and, in general, are not orthogonal. The property of being a 
basis has been proved for the eigenfunction sets of only some one-dimensional nonlinear 
problems [21,22]. However, it is important to stress that we do not need and shall not 
use these properties in what follows. 

In order to induce intermode transitions, one has to impose an additional time- 
dependent potential V = V(r,t) and to consider a coherent field ip = ip(r,t) satisfying 
the temporal nonlinear Schrodinger equation 

ih ^ = [H(<p) + v] v . (3) 

Supposing that at the initial time all atoms were condensed to the ground-state coher- 
ent mode, one has the initial condition 

V?(r,0) = v?o(r). (4) 

For transfering atoms from the ground state to a higher mode requires to impose an 
alternating potential 

V — Vi(r) cosut + V^r) sino;£ , (5) 

with a frequency uo being in resonance with the chosen transition frequency. If the 
wanted excited mode has the energy E k , the transition frequency is 

co k = j. (E k - E ) . (6) 
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Then the resonance condition tells that the detuning of u from u k is to be small, 



< 1 (AiJ = LO-LO k ) . (7) 



One expects that, under the resonance condition (7), only the considered modes, 
connected by the resonance frequency (6), will be mainly involved in the process of 
excitation. This can be proved explicitly by invoking the method of averaging [23]. 
For this purpose, let us look for the solution of equation (3) in the form 

(p(r, t) = 2 c n (t)ip n (r) exp ^ E n tj , (8) 



where c n (t) is a slowly varying amplitude, such that 

< 1 . (9) 



E n 



dc n 



dt 



The presentation (8) is to be substituted into equation (3), which is multiplied by 
exp(j-E n t) and whose right-hand side is averaged over time. Integrating over time, 
the amplitudes c n , according to condition (9), are treated as quasi-invariants. The 
integration of exponents yields 



lim ^ f exp (E m - E n )t \ dt = S mn , 
n ^oo i Jo In J 

lim ^ / exp (E m + E k - E n - E t )t } dt = 5 mn 5 k i + 5 m i5 nk - 5 mk 5 kn 
n -><x> 1 Jo In J 



T 

I r' 

im 

T 

This procedure results in the equation 



^ = -i J2 a nm\c m \ 2 c n - l - 5 n0 f3 0k c k e tAuJt - J 5 n kPo k coe lAuJt , (10) 

dt TT n 2 2 

in which the notation is used for the nonlinear transition amplitude 

N f 

a nm = A— J |v?„(r)| 2 (2|v? m (r)| 2 - |y?„(r)| 2 ) rfr , (11) 

due to atomic interactions, and for the linear amplitude 

Pok = U <P* (r) [Ki(r) - iV 2 {v)\ <p k (r) dr , (12) 

caused by the resonant field (5). Note that the orthogonality of the modes <p n (r) has 
nowhere been required. The initial condition to equation (10), in line with condition 
(4), is 

c„(0) = 5 n0 . (13) 
In the case when n ^ 0, k, the solution to equation (10) can be written as 



c n (t) = c n (0) exp < -% a nm f |c m (t')| 2 dt' 
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This, in compliance with the initial condition (13), shows that c n (t) = for all n 7^ 0, k. 
So that only the behaviour of c (t) and Ck(t) is nontrivial, with the related initial 
conditions 

c (0) = 1 , c fc (0) = . (14) 
Thus, equation (10) reduces to the system of equations 

d°0 I 1 2 in iAcut 

— = -ta ok \c k \ c - - PofcC fc e 

*^ Cfc • I 1 2 * ,o* -iAuit / 1r \ 

— = -2a fc0 |c | c fe - - /%c e , (15) 

with the initial conditions (14). The solutions to equations (15) are called guiding 
centers. 



3 Change of Variables 

Equations (15) can be simplified by changing the variables. To this end, it is convenient 
to introduce the population difference 

s = \ck\ 2 - |c | 2 • (16) 

The amplitudes Co and c k can be presented as 

i-s\ 1/2 r / Alo 



Co = ) ex P Y ( n ° + ~ t 

c k = (-y^) ex P {* (^i - *)} > ( 17 ) 

with 7r = 7To(i) and 7Ti = 7Ti (t) being real phases. 
Let us define the combination 

a = \ ("ofc + «feo) (18) 
of the amplitudes (11), which is a real quantity, present the amplitude (12) as 

A)* = /3e* (/3=|A) fc |), (19) 
and also define the renormalized detuning 

5 = Au + i (a ofe - a k0 ) . (20) 
Lastly, we introduce the phase variable 

X = Til - 7T + 7 . (21) 

The new functional variables (16) and (21) change in the rectangle 

- 1 < s < 1 , 0<x <2tt . (22) 
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The related initial conditions are 

s(0) = -1 , x(0) = x . (23) 

The value x = 717(0) — 7r (0) +7 can be any in the interval [0, 2n}. This is because, even 
if the initial phases of the considered modes were the same, the quantity 7 can be varied 
by choosing an appropriate alternating potential (5). Also, even when 7Ti(0) = 7r (0), 
the time dependence of iri(t) and 7r (t) is different, so that the evolution of x = x(t) is 
not trivial. 

With the new variables (16) and (21), equations (15) can be transformed to the 
Hamiltonian form characterized by the Hamiltonian 

H(s, x) = ^ s 2 - (3Vl - s 2 cos x + 5s . (24) 

The autonomous equations 

ds_ _dH_ dx _ OH 

dt dx ' dt ds 

are identical to those that follow from equations (15) after substituting there expres- 
sions (17), and which are 

ds r- - dx (3s . . . 

— = — pv 1 -s z sini, — = as H — cosx + . (25) 

dt dt vl — s 2 

Equations (25) are more convenient for analyzing than those (15). 

The autonomous Hamiltonian form of the evolution equations (25) tells us two 
things. First, there is no dissipation for the coherent modes. The absence of intrinsic 
decoherence, connected with dissipation, here is quite clear. For such a decoherence to 
arise, the Bose system is to be in contact with an external bath [24,25] whose role can be 
played, e.g., by largely detuned external laser beams [26] or by disturbing measurement 
instruments [27]. Decoherence of a coherent wave packet may appear owing to the 
number-of-particle fluctuations [28], which, actually, presupposes the existence of an 
external bath. The latter is a general cause of dissipation for statistical systems [29,30] 
(other references on dissipative systems, can be found in the recent review [31]). 

Another conclusion which results from the Hamiltonian representation of the au- 
tonomous evolution equations (25) is that in this dynamical system there can be no 
chaos. Hence the critical phenomena discovered [18,19] in the dynamics of fractional 
populations cannot be attributed to the appearance of chaotic motion. The origin of 
these critical phenomena will be elucidated in the following section. 



4 General Analysis 

To understand what happens with the solutions to the evolution equations (25) under 
varying parameters, we have to analyze the general phase structure of these equations. 
For simplifying formulas, it is useful to introduce the dimensionless parameters 

b=P, e= 5 -. (26) 
a a 
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The first of the latter describes a relative intensity of the alternating field (5), while the 
second parameter characterizes a relative value of the detuning (20). These parameters 
can be positive as well as negative. The sole thing is that the relative detuning will be 
treated as a small parameter, |e| <C 1. 

Due to the existence of the integral of motion (24), the trajectory, starting at the 
initial point (23) and defined by the equality H(s,x) = H(so,x ), writes 

— — bVl — s 2 cosx + es = - —e. (27) 
The right-hand sides of equations (25), with the notation (26), can be presented as 

Vi = — bVl — s 2 sin a; , v 2 = s H — cos a; + e . (28) 

Vl — s 2 

The stationary solutions to equations (25), given by v-y — v 2 — 0, are defined by the 
equations 

s 4 + 2es 3 -(l-b 2 -e 2 )s 2 -2es-e 2 = 0, sina; = 0. (29) 

The following analysis, keeping in mind the smallness of the detuning, will be accom- 
plished in the linear approximation with respect to e. Also, one has to always remember 
that physically admissible fixed-point solutions are only those that are in the rectangle 
(22). 

When b 2 > 1, there exist the fixed points 

p 

s 1 — — , x x — u , 

s* 2 = - - , x\ = 7T , (30) 

s 3 = s 1 , x 3 = zn . 
If < b < 1, the fixed points (30) continue to exist, but there arise two new points 

U2 



b 2 e 



Sl = Vl - b 2 + y— — , X\ = 7T , 

^ = -V / T^+ T ^, XI = 7T. (31) 

And if — 1 < b < 0, then the fixed points (30) again exist, but there appear the 
additional points 

si = v / r^+ T ^, x* A = o, 

Sg — S 4 , — Z7T , 

s* = s* 5 , 4 = 2tt . (32) 
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In this way, the value b 2 = 1 corresponds to a dynamical phase transition, when the 
structure of the phase portrait changes qualitatively. 

To analyse the stability of motion, we consider the Jacobian expansion matrix 
X = [Xij] with the elements 



bs 



dv 1 



SIM , 



X12 = 



dv 1 
dx 



—by/l — s 2 cosx , 



X21 = 



dv 2 
ds 



1 + 



b cos x 



1 



2^3/2 ' 



X 22 = 



dv 2 
dx 



-X 
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The local expansion rate [32], defined as 



A(t) = - Re /*Tr X(t') dt' , 
t Jo 

nullifies, A(t) = 0, as it should be for Hamiltonian systems whose phase volume is 
conserved. The eigenvalues of the expansion matrix X are given by the equation 



X 2 



b 2 ( 2 . 



sin 2 x — cos 2 x 



bVi- 



s z cosx . 



(33) 



Evaluating the eigenvalues at the fixed points, we shall employ expressions (30) to (32), 
limiting ourselves by the linear in e approximation. 
In the case of b 2 > 1, we have 



Xf = ±iy/b(b+l) = X± , Xf = ±iy/b(b-l) , 



(34) 



so that all fixed points (30) are centers. 

When < b < 1, the first and third fixed points remain the centers, while the 
second fixed point (s|, x* 2 ) becomes a saddle. The fixed points (31) are also the centers. 
The related exponents are 



X 1 ± = ±* V /&(1 + 6)=X 3 ± , 
(2 + b 2 )e ll/2 



X. 



± 



±Jb(l-b), 



1 + 



(1 - b 2 f' 2 



(2 + b 2 )e 

(1 - b 2 y>/ 2 



1/2 



(35) 

The saddle separatrices which are the trajectories that pass through the saddles and 
separate the phase regions with qualitatively different dynamic properties, are given 
by the equation H(s,x) = H(s\,x*£), which results in the separatrix equation 



— — bVl — s 2 cosx + es = b , 



(36) 



defining two separatrices, since this is a square equation with respect to s. The sepa- 
ratrix extremal points can be found from the equation 



ds 
dx 



6(1 



s sin x 



(s + e) y/l — s 2 + bs cos x 



= . 
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As follows from the above equations, the lower separatrix parts touch the boundary at 
the points s — — 1, x — 0, 2ir, provided that 

b + e = \. (37) 

In this way, accepting as initial conditions s = —1, x = 0, one has the following 
picture. For b < 0.5 — e, the motion is limited from above by the lower separatrix parts 
and from below, by the boundary s = —1. When the pumping parameter b > 0.5 — e, 
the trajectory passes to the phase region limited from above by the upper separatrix 
parts and from below by the lower separatrix parts. If b = 0.5 — e, the dynamical 
system is structurally unstable with respect to small variations of initial conditions. 
On the manifold of the system parameters, the line (37) plays the role of a critical line. 
In the vicinity of this line, solutions to the evolution equations display dramatic effects, 
when a tiny variation of a parameter qualitatively changes the properties of solutions. 
This makes it possible to call such dynamical effects the critical dynamic phenomena 
[33]. 

For the case when — 1 < b < 0, the expansion exponents, corresponding to the fixed 
points (30) and (32), are 



xt = ±J\b\{l + b) = X 3 ± , X± = ±iJ\b\(l - b) , 



xf = ±iVT^v 



xf = ±iVT^¥ 



i + 



l - 



b 2 )e 



(1 - b 2 fl\ 

(2 + b 2 )e 
(1 - b 2 f/ 2 



1/2 



1/2 

= Xf . (38) 



Hence, the first and third fixed points become the saddles, while all other points are 
the centers. The separatrices connecting the saddles are defined by the equations 
H(s,x) = H(sl,xl) = H (S3, Xg), which yield 



s 2 



— byl — s 2 coax + es = —b . (39) 



The lower separatrix part touches the boundary at the point s = — 1, x — ir under the 
condition 

\b\+e = ±. (40) 

The phase picture, as compared to the previous case < b < 1, looks similar, but 
being shifted by 7r along the axis x. Now, if the initial point would be Sq = —1, Xq = 0, 
the motion would be always limited from above by the lower separatrix parts and from 
below, by the boundary s = — 1. No dramatic changes would happen when varying 
b. However, if the initial point is taken as s = — 1, x = n, one again encounters the 
same critical dynamic phenomena on the critical line (40). 

This analysis explains that the occurrence of critical dynamic phenomena is caused 
by the existence of a critical line on the manifold of system parameters and happens 
only under a special choice of initial conditions, when the latter are touched by a 
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separatrix. The initial conditions for the variable (21) can be varied in a wide diapason 
by choosing the appropriate alternating field (5), which would yield the related linear 
amplitude (12), with the required value of 



7 = arg [3 0k = tan 



1 Im f3 0k 
Re (3 0k ' 



defined in equation (19). 



5 Numerical Solution 



In order to explicitly illustrate the critical dynamic phenomena, occurring when cross- 
ing the critical line on the parametric manifold, we have accomplished numerical cal- 
culations of the fractional populations 



This can be done by numerically solving either equations (15) or (25), which are equiv- 
alent. The results are presented in figures 1 to 3, where time is measured in units of 
a -1 , the pumping parameter b = 0.5 is fixed, while the detuning e is varied so that to 
cross the critical line (37). The initial conditions for the evolution equations are taken 
as s = — 1, = 0, that is n (0) = 1, n k (0) = 0. 

In figure 1, the detuning is negative, so that one is slightly below the critical line 
(37). The fractional populations oscillate displaying a kind of Rabi oscillations, if one 
looks for analogies with optics [34]. However, since the problem considered here is 
nonlinear, there is no a well defined constant Rabi frequency, whose analog would be 
now a function of time [18,19]. In some sense, these oscillations could be named non- 
linear Rabi oscillations. Approaching the critical line (37) by increasing the detuning 
e, the oscillation amplitude increases. The motion of the population n k of the excited 
topological mode is limited from above by the lower parts of separatrices and from 
below, by the boundary n k = 0. 

Figure 2 demonstrates the motion in a tiny vicinity of the critical line (37). Chang- 
ing the detuning from e = 10~ 4 to e = drastically transforms the shape of oscillations. 
The oscillation period more than doubles and wide flat platos in the time dependence 
of the population appear. Crossing the critical line, with a microscopic variation of 
the detuning from e = to e = 10~ 9 , again yields a drastic transformation of the 
population shapes. The period is again almost doubled; the upward casps of n k (t) 
and down- ward casps of n (t) arize. The appearance of these casps means that the 
motion has passed to another phase region, as has been discussed in the analysis of the 
previous section. 

After crossing the critical line (37), the dynamics of the fractional populations re- 
mains qualitatively unchanged. Increasing the detuning slightly changes the oscillation 
period and smoothes the shape of the oscillation curves, as is shown in figure 3. 

In conclusion, we have presented a theory for the resonant excitation of topological 
coherent modes of trapped Bose atoms. These modes are the stationary solutions to 



ra (t) = |c (t)| 2 



1 - s(t) 



n k {t) = \c k (t)\ 2 



1 + s(t) 
2 



2 
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the nonlinear Schrodinger equation, which is also sometimes called the Gross-Pitaevskii 
equation. This equation provides a correct description of trapped atoms at low tem- 
peratures [1-3,35]. The principally important part of the paper is the demonstration 
of the occurrence of critical dynamic phenomena in the process of exciting coherent 
modes and a thorough elucidation of the origin of these phenomena. 
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Figure Captions 



Figure 1. The fractional populations of the ground coherent mode no(t) (dashed 
line) and of the excited topological coherent mode n k (t) (solid line) as functions of 
time, measured in units of a' 1 , for the fixed pumping parameter b = 0.5 and the 
negative detuning parameter below the critical line: (a) e = — 1CT 1 ; (b) e = — 1CT 2 ; (c) 

£= -10" 4 . 



Figure 2. Dramatic changes in the dynamics of the fractional populations of the 
ground coherent mode (dashed line) and excited mode (solid line) when crossing the 
critical line on the parametric manifold for fixed b = 0.5 and varying detuning: (a) 
e = 0; (b) e = 10~ 9 . 

Figure 3. Qualitatively different temporal behaviour of the fractional populations 
of the ground mode (dashed line) and excited mode (solid line) after crossing the 
parametric critical line, with fixed b = 0.5 and varying detuning: (a) e = 10 -4 ; (b) 
£ = to- 2 ; ( c ) e = 10" 1 . 
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